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Abstract 

We show that within a recently developed nonlocal chiral quark model the critical density for a 
phase transition to color superconducting quark matter under neutron star conditions can be low 
enough for these phases to occur in compact star configurations with masses below 1.3 M@. We 
study the cooling of these objects in isolation for different values of the gravitational mass. Our 
equation of state (EoS) allows for 2SC quark matter with a large quark gap ~ 100 MeV for u and 
d quarks of two colors that coexists with normal quark matter within a mixed phase in the hybrid 
star interior. We argue that, if the phases with unpaired quarks were allowed, the corresponding 
hybrid stars would cool too fast. If they occured for M < 1.3 Mq, as it follows from our EoS, 
one could not appropriately describe the neutron star cooling data existing by today. We discuss a 
"2SC+X" phase, as a possibility to have all quarks paired in two-flavor quark matter under neutron 
star constraints, where the X-gap is of the order of 10 keV - 1 MeV. Density independent gaps do 
not allow to fit the cooling data. Only the presence of an X-gap that decreases with increase of the 
density could allow to appropriately fit the data in a similar compact star mass interval to that 
following from a purely hadronic model. This scenario is suggested as an alternative explanation 
of the cooling data in the framework of a hybrid star model. 



I. INTRODUCTION 



Recently, we have reinvestigated the cooling of neutron stars within a purely hadronic 
model Q], i.e., suppressing the possibility of quark cores in neutron star interiors. We have 
demonstrated that the neutron star cooling data available by today can be well explained 
within the "Nuclear medium cooling scenario", cf. , i.e., if one includes medium effects 
into consideration. This scenario puts a tight constraint to the density dependence of the 
asymmetry energy in the nuclear equation of state (EoS). The latter dependence is an 
important issue for the analysis of heavy ion collisions especially within the new CBM 
(compressed baryon matter) program to be realized at the future accelerator facility FAIR 
at GSI Darmstadt. The density dependence of the asymmetry energy determines the proton 
fraction in neutron star matter and thus governs the onset of the very efficient direct Urea 
(DU) process. This process, once occurring, would lead to a very fast cooling of neutron stars. 
Thereby, if neutron stars with M < cooled slowly, the stars with the mass M only 

slightly above M^Y would be already characterized by a very rapid cooling, yielding surface 
temperatures significantly below those given by the X-ray data points. Unfortunately, the 
mass spectrum of cooling neutron stars is not known. One can not rely on neutron star 
mass measurements in binary radio pulsars (BRPs) with the averaged neutron star mass 
value M B rp = 1.35 ± 0.04 M & J, [| because these objects come from a twice selected 
population: probably not all neutron stars go through the active radio-pulsar stage p , and 
the properties of neutron stars in binaries may be different from those of isolated objects, 
see e.g. jzj]- However, it seems not a fair assumption that most of cooling neutron stars 
should be essentially lighter than the lower limit for BRP masses MBRp, m i n = 1-31 M & . 
Therefore, the threshold for the DU process should better occur at high enough densities 
such that typical neutron stars are not affected by it. This assumption about the mass 
distribution can be developed into a more quantitative test of cooling scenarios when these 
are combined with population synthesis models. The latter allow to obtain Log N - Log S 
distributions for nearby coolers which can be tested with data from the ROSAT catalogue 
^]. In addition the microscopically based Argonne model for the EoS, that is supposed to 
be the most realistic one j^, yields a critical mass for the DU process ~ 2 M & . At these high 
star masses the central baryon density exceeds already the 5-fold of the nuclear saturation 
density such that exotic states of matter as, e.g., hyperonic matter or quark matter are more 
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appropriate. Ref. argued that the presence of quark matter in massive compact star 
cores is a most reliable hypothesis since within their model hyperons appear at sufficiently 
low densities and the equation of state becomes too soft to carry neutron star masses in the 
range of M BRP Q. 

For quark matter in neutron stars, the interaction (a s ^ 0) permits the quark DU process. 
The works 3, Q Q have shown that this result does not exclude the possibility that 
neutron stars might possess large quark matter cores. The suppression of the emissivity of 
the DU process may arise due to the occurence of diquark pairing and thus color supercon- 
ductivity in low-temperature quark matter. To demonstrate this possibility one has used the 
simplifying assumption that all quarks are gapped with the same gap which is not valid for 
superconducting phases at low densities. In the present paper we abandon that restriction 
and suggest a more subtle pairing pattern that is constrained by comparison with neutron 
star cooling phenomenology. 

II. COLOR SUPERCONDUCTIVITY 

The quark-quark interaction in the color anti-triplet channel is attractive driving the 
pairing with a large zero-temperature pairing gap A ~ 100 MeV for the quark chemical 



potential \i q ~ (300 -r- 500) MeV, cf. [l6|, Il7| . for a review see |3] and references therein. 



The attraction comes either from the one-gluon exchange, or from a non-perturbative 4- 
)oint interaction motivated by instantons 19], or from non-perturbative gluon propagators 
2p| . Various phases are possible. The so called 2-flavor color superconductivity (2SC) 

phase allows for unpaired quarks of one color, say blue. There may also exist a color-flavor 



locked (CFL) phase 21| for not too large values of the dymanical strange quark mass or 



in other words for very large values of the baryon chemical potential |22j, where the color 
superconductivity is complete in the sense that the diquark condensation produces a gap for 
quarks of all three colors and flavors. The value of the gap is of the same order of magnitude 
as that in the two-flavor case. There exist other attractive quark pairing channels. However 
the attraction is rather weak. One of these channels allows for a color spin locking (CSL), 
where all quarks are paired in a spin-1 channel. The weak pairing channels are characterized 
by gaps typically in the interval ~ 10 keV -f-1 MeV. 
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III. HYBRID STARS 

In describing the hadronic part of the 
V18 + 5v + UIX* model of the EoS given in 



ivbrid star, as in 111. we exploit the Argonne 
which is based on the most recent models for 



the nucleon-nucleon interaction with the inclusion of a parameterized three-body force and 
relativistic boost corrections. Actually we continue to adopt an analytic parameterization of 
this model by Heiselberg and Hjorth- Jensen 



23| . hereafter HHJ, that smoothly incorporates 
causality at high densities. The HHJ EoS fits the symmetry energy to the original Argonne 
V18 + Sv + UIX* model in the mentioned density interval yielding the threshold density for 
the DU process n? u ~ 5.19 n (M DU ~ 1.839 M ra ). 

nn 

The 2SC phase occurs at lower baryon densities than the CFL phase, see |24J, 1251 ] . For 
applications to compact stars the omission of the strange quark flavor is justified by the fact 
that chemical potentials in central parts of the stars do barely reach the threshold value at 
which the mass gap for strange quarks breaks down and they appear in the system 



271 ] for the case 



appear a mixed phase region, see 
tension and charge screening. Refs 



We will employ the EoS of a nonlocal chiral quark model developed in 
of neutron star constraints with a 2SC phase. It has been shown in that work that the 
Gaussian formfactor ansatz leads to an early onset of the deconfinement transition and such 
a model is therefore suitable to discuss hybrid stars with large quark matter cores [3] . 
In some density interval above the onset of the first order phase transition, there may 

29!. Ref. disregarded finite size effects, such as surface 
] on the example of the hadron-quark mixed phase have 
demonstrated that finite size effects might play a crucial role substantially narrowing the 
region of the mixed phase or even forbidding its appearance in a hybrid star configuration. 
Therefore we omit the possibility of the hadron-quark mixed phase in our model assuming 
that the quark phase arises by the Maxwell construction. In this case we found a tiny density 
jump on the boundary from 7i|? adr ~ 0.44 fm~ 3 to 7?q uark ~ 0.46 fm~ 3 . 

A large difference between chemical potentials of u and d quarks forbids the pure 2SC 
phase, cf. 27]. The CFL phase is still not permitted at such densities. Ignoring weak 
pairing patterns there could be two possibilities: either the quark matter is in the normal 
phase, or there appears a region of the 2SC - normal quark mixed phase. Disregarding finite 
size effects, as in (29 1, Refs. Q, |3] have argued for a broad region of the 2SC - normal 
quark mixed phase instead of a pure 2SC phase. In our model the mixed phase continues 
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up to the central density n cen t T . Thus we have the mixed phase in the density interval: 
n hadr < n^ aix ' 1 ^ < ^quark < n < n ^"- %x ^ — n centr . In the given case the arguments of [3^ 
might be partially relaxed since the surface tension for the 2SC - normal quark boundary 
should be significantly smaller than for the quark - hadron boundary. Indeed, as has been 
indicated in Ref. 15[, the surface tension is proportional to the difference of the ener gies 
of two phases, being proportional to (A//z g ) 2 C 1 in the given case. Recently, Ref. [31| 
estimated the surface tension to be typically smaller than 5 MeV/fm 2 for densities of our 
interest. For the choice of the Gaussian formfactor the mixed phase appears for sufficiently 
low density n > n { ™ ix ' 1] ^ ^ adr = 0.44 fm -3 (M > M c quark = 1.214 M Q ). The presence or 
absence of the 2SC - normal quark mixed phase instead of only one of those phases is not 
so important for the hybrid star cooling problem since the latter is governed by processes 
involving either normal excitations or excitations with the smallest gap. 

In Fig. E]we present the mass-radius relation for hybrid stars with HHJ EoS vs. Gaussian 
nonlocal chiral quark separable model (SM) EoS. Configurations with the 2SC - normal quark 
mixed phase, "HHJ-SM with 2SC" (HHJ hadronic EoS together with Gaussian nonlocal 
chiral quark model SM permitting 2SC - normal quark mixed phase) given by the solid line, 
are stable, whereas without 2SC large gap color superconducting matter ("HHJ-SM without 
2SC") no stable hybrid star configuration is possible. In the case "HHJ-SM with 2SC" the 
maximum neutron star mass proves to be 1.793 M Q . For an illustration of the constraints 
on the mas-radius relation which can be derived from compact star observations we show 
in Fig. ^the compactness limits from the thermal emission of the isolated neutron star RX 

n 

J1856. 5-3754 as given in [22j and from the redshifted absorption lines in the X-ray burst 
spectra of EXO 0748-676 given in [33]. These are, however, rather weak constraints. 

Additionally, within the "HHJ-SM with 2SC" mixed phase we will allow for the possibility 
of a weak pairing channel for all the quarks which were unpaired, with typical gaps Ax ~ 
10 keV -f-1 MeV, as in the case of the CSL pairing channel, see jsJ, Q Since we don't 
know yet the exact pairing pattern for this case, we call this hypothetical phase "2SC+X" . 
In such a way all the quarks get paired, some strongly in the 2SC channel and some weakly 
in the X channel. 
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IV. COOLING 



For the calculation of the cooling of the hadron part of the hybrid star we use the same 
model as in The main processes are the medium modified Urea (MMU) and the pair 
breaking and formation (PBF) processes. The HHJ EoS is adopted. Below we choose two 
cooling models used in jjj distinguished by different sets of the nucleon gaps. The first choice 
(further model I) is shown in Fig. 2 (Fig. 12 of Ref. ujT). We use a fit for the relation between 
surface and crust temperatures (called "our fit" in l|). We also made calculations for the 
other two crust models used in Ref. Q]. Although the results are affected by the choice of 
the crust model, our conclusions remain untouched. The possibilities of pion condensation 
and of other so called exotic processes are suppressed since in the model 1| these processes 
may occur only for neutron star masses exceeding M^ uark = 1.214 M & . The DU process is 

irrelevant in this model up to very large neutron star mass M > 1.839 M fl . The lSn neutron 

o o 

and proton gaps are taken the same as those shown by thick lines in Fig. 5 of Ref. [1[. We 
pay particular attention to the fact that the 3P2 neutron gap is additionally suppressed by 
the factor 0.1 compared to that shown in Fig. 5 of jl}. The latter suppresion is motivated 
by the result of the recent work 3(| and is required to fit the cooling data. 

For the calculation of the cooling of the quark core in the hybrid star we use the model 
Q|. We incorporate the most efficient processes: the quark direct Urea (QDU) processes 
on unpaired quarks, the quark modified Urea (QMU), the quark bremsstrahlung (QB), the 
electron bremsstrahlung (EB), and the massive gluon-photon decay (see [3]). Following 
we include the emissivity of the quark pair formation and breaking (QPFB) processes. The 
specific heat incorporates the quark contribution, the electron contribution and the massles 
and massive gluon-photon contributions. The heat conductivity contains quark, electron 
and gluon terms. 

We are basing on the picture presented in Fig. El and add the contribution of the quark 
core. For the Gaussian formfactor the quark core occurs already for M > 1.214 M & according 
to the model j^, see Fig- d Most of the relevant neutron star configurations (see Fig. EJ) 
are then affected by the presence of the quark core. First we check the possibility of the 
2SC+ normal quark phases. 

Fig. El shows the cooling curves calculated with the given Gaussian ansatz. The variation 
of the gaps for the strong pairing of quarks within the 2SC phase and the gluon-photon 
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mass in the interval A,m s _ 7 ~ 20 -j- 200 MeV only slightly affects the results. The main 
cooling process is the QDU process on normal quarks. We see that the presence of normal 
quarks entails too fast cooling. The data could be explained only if all the masses lie 
in a very narrow interval (1.21 < M/M & < 1.22 in our case). In case of the other two 
crust models the resulting picture is similar. The existence of only a very narrow mass 
interval in which the data can be fitted seems us unrealistic as by itself as from the point 
of view of the observation of the neutron stars in binary systems with different masses, e.g., 
M B i9i3+i6 ^ 1-4408 ± 0.0003 M© and Mj 737_3039b ^ 1-250 ± 0.005 M Q , cf. Thus the 
data can't be satisfactorily explained. 

In Fig. 0]we permit the weak pairing pattern for all quarks which before in Fig. 3 were 
assumed to be unpaired. We use Ax — 1 MeV for the corresponding quark gap. Fig. 
H] demonstrates too slow cooling. Even slow cooling data are not explained. The same 
statement holds for all three crust models. Thus the gaps for formerly unpaired quarks 
should be still smaller in order to obtain a satisfactory description of the cooling data. 

In Fig. El we again allow for the weak pairing pattern of those quarks which were assumed 
to be unpaired in Fig. El but now we use Ax = 30 keV for the corresponding residual 
quark pairing gap (X). We see that all data points, except the Vela, CTA 1 and Geminga, 
correspond to hybrid stars with masses in the narrow interval M = 1.21 -r- 1.22 M , very 
similar to that in Fig. El This also seems to be unrealistic with the same argument as for the 
presence of unpaired quarks, see Fig. 3 above. We have obtained a similar picture for the 
other two crust models. We have varied the constant A in wide limits but the qualitative 
picture does not change. 

Therefore we would like to explore whether a density-dependent X-gap could allow a 
description of the cooling data within a larger interval of compact star masses. We don't 
choose the X-gap as a rising function of the chemical potential since the feature of such 
a model will not spread the accessible mass interval satisfying the cooling data constraint 
relative to the density independent X-gap case. We employ the ansatz 

Ax(jti) = A c exp[-a(/i - /i c )/// c ] , (1) 

where the parameters are chosen such that at the critical quark chemical potential fi c = 330 
MeV for the onset of the deconfinement phase transition the X-gap has its maximal value 
of A c = 1.0 MeV and at the highest attainable chemical potential /i max = 507 MeV, i.e. in 
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the center of the maximum mass hybrid star configuration it falls to a value of the order of 
10 keV. We choose the value a = 10 for which Ax(^ mas ) = 4.6 keV. In Fig. E]we show the 
resulting cooling curves for the ansatz We observe that the mass interval for compact 
stars which obey the cooling data constraint ranges now from M = 1.32 M & for slow coolers 
up to M = 1.75 Mq for fast coolers such as Vela, cf. with that we have found with the 
purely hadronic model [lj with different parameter choices. Note that according to a recently 
suggested independent test of cooling models by comparing results of a corresponding 
population synthesis model with the Log N - Log S distribution of nearby isolated X-ray 
sources the cooling model I did not pass the test. Thereby it would be interesting to see 
whether our quark model within the gap ansatz ([!} may pass the Log N - Log S test. 

In different hadronic cooling models have been presented which can explain the avail- 
able cooling data by varying the neutron star mass in some limits. The question arises 
whether one of these alternatives could be preferred. In Ref. Q we conjectured that the 
model based on Fig. 20 of [lj (Fig. of the given work) that has used another choice of 
the gaps demonstrates the best fit since it shows a more regular dependence on the neutron 
star masses. We call it model II. In this model II the ISq neutron gap is the same as in 



the model I used above, whereas the proton ISq gap is calculated by 



(see thin lines in 



381 ] but then is additionally suppressed 



Fig. 5 of Q). The neutron 3P2 gap is taken from 
by factor 0.1 as required by the result of and in order to appropriately fit the cooling 
data. According to [8] this model II succesfully passed the mentioned Log N - Log S test 

n n 

supporting the conjecture of |8j and an earlier conjecture of |3j| that neutron star masses 
should be essentially different. 

We have performed the same hybrid cooling calculation as described above, now for the 
model II. The results are qualitatively the same as those presented in Figs H3 El The 
difference is that now the curves related to low mass objects (M < M™) are pushed 
downwards describing the Crab data point. However, as in model I, both scenarios either 
that there are unpaired quarks and that the X-gaps are density independent, do not allow 
to appropriately fit the data. 

The results for the density-dependent X-gap according to Eq. (Q) now using the hadronic 
cooling model II (demonstrated by Fig. [7J) are shown in Fig. |HJ In comparison with the 
result of Fig. E] the variation of the cooling data between slow and fast coolers is obtained 
by choosing hybrid star masses in a wider range 1.1 < M/M & < 1.7. Next, it would be 
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interesting to check whether this model passes the Log N - Log S test. 
V. CONCLUSION 

Concluding, we demonstrated that within a recently developed nonlocal, chiral quark 
model the critical densities for a phase transition to color superconducting quark matter 
can be low enough for these phases to occur in compact star configurations with masses 
below 1.3 Mq. For the choice of the Gaussian formfactor the 2SC - normal quark matter 
mixed phase arises at M ~ 1.21 M Q . We have shown that without a residual pairing the 
2SC quark matter phase could describe the cooling data only if compact stars had masses 
in a very narrow band around the critical mass for which the quark core can occur. Since 
there are observations of neutron stars with higher and essentially different masses such a 
scenario should be disfavored. 

Then we assumed that formally unpaired quarks can be paired with small gaps Ax < 
1 MeV (2SC+X pairing), which values we varied in wide limits. With density independent 
gaps we failed to appropriately fit the cooling data. 

We showed that the present day cooling data could be still explained by hybrid stars, 
however, when assuming a complex pairing pattern, where quarks are partly strongly paired 
within the 2SC channel, and partly weakly paired with gaps Ax < 1 MeV, being rapidly 
decreasing with the increase of the density. Within the given hypothesis and in the framefork 
of our quark model for the equation of state the observational data are explained in a broad 
mass interval M — 1.32 - 1.75 M for the model I and M — 1.1 - 1.7 M for the model 
II, as it is favored by the experimental data. However we will point out that the result is 
model dependent and demonstrates only a possibility that the presence of the 2SC+X phase 
might not contradict to the neutron star cooling data. 

We have found rather strong constraints for the density dependent X-gap pairing. It 
remains to be investigated which microscopic pairing pattern could fulfill the constraints 
obtained in this work. Another indirect check of the model could be the Log N - Log S test. 
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FIG. 1: Mass - radius relations for compact star configurations with different EoS: purely hadronic 
star with HHJ EoS (dashed line), stable hybrid stars with HHJ - Gaussian nonlocal chiral quark 
separable model (SM) with 2SC phase (solid line) and with HHJ - SM, without 2SC phase (dash- 
dotted line). Data for two sources are also indicated, see text. 
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G. 2: Model I. Cooling curves according to the nuclear medium cooling scenario, see Fig. 12 of 
lj. The labels correspond to the gravitational masses of the configurations in units of the solar 



mass. 
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FIG. 3: Model I. Cooling curves for hybrid star configurations with Gaussian quark matter core 
in the 2SC phase. The labels correspond to the gravitational masses of the configurations in units 
of the solar mass. 
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FIG. 4: Model I. Cooling curves for hybrid star configurations with Gaussian quark matter core 
in the 2SC+X phase. The weak pairing gap is Ax — 1 MeV. The labels correspond to the 
gravitational masses of the configurations in units of the solar mass. 
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FIG. 5: Same as Fig. 0]but with a weak pairing gap of 30 keV. 
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FIG. 6: Same as Fig. 0]but with a density dependent pairing gap according to Eq. ([T|). 
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FIG. 7: Model II. Same as Fig. 20 of Ref. 



D. 
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FIG. 8: Same as Fig. El but for model II. 
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